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Abstract 



00 , 

, We analyze the quantum dynamics of radiation propagating in a single mode optical 

fiber with dispersion, nonlinearity, and Raman coupling to thermal phonons. We 
start from a fundamental Hamiltonian that includes the principal known nonlinear 
effects and quantum noise sources, including linear gain and loss. Both Markovian 
and frequency-dependent, non-Markovian reservoirs are treated. This allows quan- 
CN ! turn Langevin equations to be calculated, which have a classical form except for 

additional quantum noise terms. In practical calculations, it is more useful to trans- 
form to Wigner or +P quasi-probability operator representations. These result in 
stochastic equations that can be analyzed using perturbation theory or exact nu- 
S^h. merical techniques. The results have applications to fiber optics communications, 

networking, and sensor technology. 

a 

3. 060.4510, 270.5530, 270.3430, 190.4370, 190.5650, 060.2400 



1. Introduction 

The propagation of electromagnetic radiation through optical fibers is the central paradigm 
of optical communications and sensor technology. It is also a novel physical system, due to 
the materials processing of fused silica, that leads to single-mode behaviour with extremely 
low losses. Over short distances (depending on the pulse intensity) the well-known nonlinear 
Schrodinger equation can describe most optical fibers with great accuracy, and leads to 
soliton behaviour, as well as to many other effects. Over longer distances, a number of 
reservoir effects intervene, including attenuation, Raman scattering, and the use of amplifiers 
and filters to compensate for losses. At the quantum level, both the original nonlinearity 
and the additional couplings to reservoirs can lead to quantum noise - which modifies the 
predictions of the classical nonlinear Schrodinger equation. 

In this paper, we analyze the effects of quantum noise in fiber optics. This extends 
and explains in more detail earlier theoretical work in this area, which led to the first 
predictionEH and measurement! of intrinsic quantum noise effects in optical solitons. The 
theory given here includes a detailed derivation of the relevant quantum Hamiltonian. We 
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include quantum noise effects due to nonlinearities, Raman reservoirs and Brillouin scat- 
tering. The Raman/ Brillouin noise is modeled using a multiple Lorentzian fit to measured 
fluorescence data, in order to estimate the Raman gain coefficients. Both gain and loss 
effects are included. This treatment is unified with theory of gain/loss reservoirs, which 
was also predicted^ and observedi to have large effects on soliton propagation. All these 
reservoirs are treated without using the Markovian approximation, in order to accurately 
treat the frequency- dependent reservoirs found in practical applications. 

The purpose of this work is to lay the foundations of practical methods for calculating 
and numerically simulating quantum effects in nonlinear optical fibers. These are significant 
whenever quantum-limited behaviour is important in communications, sensing, or measure- 
ment with optical fibers. 

We introduce the basic quantum Hamiltonian for an optical fiber in section |2|. This gives 
a Heisenberg equation of motion which reduces to the nonlinear Schrodinger equation in 
the classical limit. The equation of motion is extended in section |3] to include Raman and 
Brillouin effects, with gain and absorption processes considered in section f|. The complete 
Heisenberg equation in section [5] is the central result of this paper. Stochastic partial differ- 
ential equations can be derived from the quantum equation, using the phase-space methods 
outlined in section ^. Applications of these methods to practical examples are reserved for 
a following paper (QNII)Q. 



2. Nonlinear Schrodinger Model 

The interaction between photons in a fiber is mediated through the dielectric material con- 
stituting the fiber. The coupling to the dielectric introduces frequency dependent and time 
delayed behaviour. The complete Hamiltonian and its derivation have been given in the 
literatureSEHlil; we will only go over the salient points here. The starting point is a La- 
grangian that generates the classical Maxwell's equations for a one- dimensional dielectric 
waveguide, and that gives a Hamiltonian corresponding to the dielectric energyB: 



H D = fdV — |B| 2 + /*E(t') -D(*'W 
J 2/i Jt 



(2.1) 



where the electric field E = (D — P)/eo includes the polarization response of the dielectric 
to an incident electric displacement D. The field variables are then quantized by introduc- 
ing equal-time commutators between the canonical coordinates D and B. We note that, of 
course, it is also possible to make other choices of canonical momenta. This choice corre- 
sponds to a dipole-coupled", rather than minimal-coupled fundamental Lagrangian. While 
different Lagrangians are canonically equivalent, the present choice - originally introduced^ 
by Hillery and Mlodinow in applications to dielectric theory - has the advantage of com- 
parative simplicity. The Lagrangian must produce both the correct energyEl and Maxwell's 
equations, otherwise the conjugate momenta will contain an arbitrary scaling, leading to 
incorrect commutation relationsllQ. 
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A. Fiber-optic Hamiltonian 



The optical fiber treated will be a single transverse mode fiber with dispersion and nonlin- 
earity. Since boundary effects are usually negligible in experiments, it is useful to first take 
the infinite volume limit, which effectively replaces a summation over wave- vectors with the 
corresponding integral. We will start with a single polarization direction (i.e., a polarization 
preserving fiber). The more general case is summarized elsewhere^, and will be treated in 
detail subsequently. The basic normally ordered, nonlinear Hamiltonian for the fiber in this 
case isl: 



dkhu(k)a) (k)a(k) - / d 3 x 



Ax 



2e(w ) 
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(2.2) 



Here u(k) is the angular frequency of modes with wave- vector k, describing the linear po- 
lariton excitations in the fiber, including dispersion. We will assume that u)(k) describes the 
average linear response of the fiber, in the limit of a spatially uniform environment. If the 
fiber is spatially nonuniform, then it is necessary to add additional inhomogeneous terms to 
the Hamiltonian, of generic form Ax^^(x). As usual, e(uo) is the mode-average dielectric 
permittivity at a carrier frequency ujq = w(k ), while d(k) is an annihilation operator defined 
so that 



a(fc'),a t (fc)l = 5(k - k') 



(2.3) 



The coefficient x®(x) is the nonlinear coefficient arising when the electronic polariza- 
tion field is expanded as a function of the electric displacement, in the commonly used 
Bloembergenll! notation (the units are S.I. units, following current standard usage). This 
may vary along the longitudinal position on the fiber, if the fiber has a variable compo- 
sition. In terms of modes of the waveguide, and neglecting modal dispersion, the electric 
displacement field operator D(x) is: 



D( 



dk 



hke(cu(k))v(k) 
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d(k)u(r) exp(ikx) + h.c. 



(2.4) 



(2.5) 



Here v(k) = du(k)/dk is the group velocity. The function u(r) gives the transverse mode 
structure. Although a general mode structure can be included, for the purposes of this paper 



where: 
'd 2 r|u(r)| 2 = 1 . 



we could equally well assume a square wave-guide of area A, which gives u(r) ~ e y /y A. 
We note here that the above mode expansion for a dispersive medium is a rather general 
one, and has been worked out both from macroscopic quantizationQ, and from a microscopic 
mo delS with an arbitrary number of electronic or phonon resonances. 

In the infinite volume limit, the polariton field is defined by noting that the annihilation 
and creation operators can be related to a quantum field using: 
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$(t,x) 



Z7T 



d/c a(i, fc) exp[i(k — k )x + iuot] 



(2.6) 



This photon-density operator ty(t,x) is the slowly varying field annihilation operator for 
the linear quasi-particle excitations of the coupled electromagnetic and polarization fields 
traveling through the fiber!. The nonzero equal-time commutations relations for these Bose 
operators are 



$(t,z),$t( fj 



X 



S(x 



X 



(2.7) 



As shown in earlier treatments!, the Hamiltonian [Eq. (|2.2|) 1 can now be rewritten ap- 
proximately as: 



H F = hJdxJ dx'w(x,x')&(t,x)$(t,x') - ~J dxx E {x)$ f2 {t,x)$ 2 (t,x) . 



(2.8) 



Here we have introduced the kernel u(x,x'), which is the linear dielectric component of the 
Hamiltonian, and a nonlinear coupling term xe{x)- This kernel is then Taylor expanded 
around k = k Q , and approximated to quadratic order in {k — ko), by: 



U{X, X 



J ^-u(k) exp^k - k )(x - x')} - ~k v(k Q ) J d 2 rA X «(x)|u(r)| 2 5(x 



x 



dk 



u + Au(x)]6(x -x')+ / —[iuj' (d x/ - d x ) + u'^d^) H }exp[ik(x - x')} . (2.9) 



47T 

In writing down Eq. ( |2.8|) , we have assumed that the frequency dependence in the non- 
linear coupling can be neglected, which is a good approximation for relatively narrow band- 
widths. The nonlinear term is often called the x^ effect, so named because it arises from 
the third order term in the expansion of the polarization field in terms of the electric fieldlll. 
It causes an electronic contribution n 2e to the intensity dependent refractive index, where: 
n = n + In 2 = n + I{n 2e + n 2p ). Thus we define x E > m units of [m/s], as: 



3hwlv(k ) 
4e(w )c 2 



J d 2 vx 



(3) ^x)|u(r)| 4 



hn 2e (x)u 2 v 2 
17 



(2.10) 



Here A — [J <i 2 r|u(r)| 4 ] -1 is the effective modal cross-section, and n 2e is the refractive index 
change per unit field intensity due to electronic transitions. This is less than the total value 
observed for n 2 , since phonon contributions have yet to be included. 

The free evolution part of the total Hamiltonian, which will be removed in subsequent 
calculations, just describes the carrier frequency ujq. This is not needed in Heisenberg picture 
calculations for ty(t, x), since it is cancelled by the slowly varying field definition. Next, on 
partial integration of the derivative terms and Fourier transforming, the resulting interaction 
Hamiltonian H' F describing the evolution of $ in the slowly varying envelope and rotating- 
wave approximations is: 

H' F = H F — J dkhuod) (k)a(k) 

= - dx Aw(i)f f f + — ( - ^ f W) + — W f W 

2 J —oo 2 ^ ^2 



(2.11) 
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For simplicity, only quadratic dispersion is included here. However, the extension to 
higher-order dispersion is relatively straightforward. This can be achieved by including 
higher-order terms in the expansion of the dielectric kernel, or else by treating the dis- 
persion as part of the reservoir response function - as in following sections. The response 
function approach has the advantage that a completely arbitrary polarization response can 
be included, and transformations to a different frame of reference are simplified. If part of 
the dielectric response is treated using response functions, then this part of the measured 
refractive index must be excluded from the free Hamiltonian, to avoid double-counting. 



B. Heisenberg equation 



From the interaction Hamiltonian [Eq. |2.11 |, we find the following Heisenberg equation of 



motion for the field operator propagating in the +x direction: 

§{t,x) . (2.12) 



iuj" d 2 - 
-iAw(x) + + iX (x)W(t, x)^(t, x) 



This is the quantum nonlinear Schrodinger equation in the laboratory frame of reference, 
which is completely equivalent to the theory of a Bose gas of massive quasi-particles with an 
effective mass of %/uj" and an average velocity of v , for photons near to the carrier frequency 
of interest. It includes the possibility that the dielectric constant (i.e., the linear response) 
has a spatial variation, through the term Au(x) . 

We note here that it is occasionally assumed that operators obey equal-space, rather 
than equal-time commutation relations. This cannot be exactly true, since commutators in 
an interacting quantum field theory are only well-defined at equal times. At different times, 
it is possible for a causal effect to propagate to a different spatial location, which can there- 
fore change the unequal-time commutator. In other words, imposing free-field equal-space 
but unequal-time commutators is not strictly compatible with causality. The assumption 
of equal-space commutators may be used as an approximation under some circumstances, 
provided interactions are weak. In this paper, we will use standard equal-time commutators. 



3. Raman Hamiltonian 

To the Hamiltonian given in Eq. ( p.llj ) must be added couplings to linear gain, absorption 
and phonon reservoirsEHll. The gain and absorption reservoirs are discussed at length in 
section |j. The phonon field consists of thermal and spontaneous excitations in the dis- 
placement of atoms from their mean locations in the dielectric lattice. Although previous 
quantum treatments of Raman scattering have been givenEl, it is necessary to modify these 
somewhat in the present situation. The Raman interaction energySll of a fiber, in terms 
of atomic displacements from their mean lattice positions, is known to have the form: 

//// = ^v*to&M*)M + ^X>i ■ mm • (3- 1 ) 

3 ij 

Here D(x J ) is the electric displacement at the j-th mean atomic location x J , 5x J is the atomic 
displacement operator, r)f is a Raman coupling tensor, and Kij represents the short-range 
atom-atom interactions. 
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In order to quantize this interaction with atomic positions using our macroscopic quan- 
tization method, we must now take into account the existence of a corresponding set of 
phonon operators, b{ui,x) and b\u,x). These operators diagonalize the atomic displace- 
ment Hamiltonian in each fiber segment, and have well-defined eigen-frequencies. There 
are calculations^ of the frequency spectrum and normal modes of vibration for vitreous 
silica, using physical models based on the random network theory of disordered systems. 
The computed vibrational frequency spectrum is remarkably similar to the observed Ra- 
man gain profiled. The phonon-photon coupling induces Raman transitions and scatter- 
ing from acoustic waves (the Brillouin effect) resulting in extra noise sources and an ad- 
ditional contribution to the nonlinearity. The initial state of phonons is thermal, with 
n th (u) = [exp (hu/kT) - l] -1 . 



A. Hamiltonian and Heisenberg equations 



In terms of these phonon operators, the fiber Hamiltonian in the interaction picture and 



within the rotating wave approximation for a single polarization i M H' 
we have introduced a Raman interaction Hamiltonian: 



Hr + H' F , where 



H, 



h 



dx 



dw {&(x)$(x)R(u},x) b(u,x) + V(to,x) + ujV(uj, x)b(u, x)\ . (3.2) 



Here, the atomic vibrations within the silica structure of the fiber are modeled as a continuum 
of localized oscillators, and are coupled to the radiation modes by a Raman transition with 
a real frequency dependent coupling R{lo,x). This coupling could be nonuniform in space, 
and is determined empirically through measurements of the Raman gain spectrum0. The 
atomic displacement is proportional to b + b\ where the phonon annihilation and creation 
operators, b and , have the equal-time commutations relations 



b(t, u, x),tf(t, lo', x') = 5(x — x')S(u — to') 



(3.3) 



Thus the Raman excitations are treated as an inhomogeneously broadened continuum of 
modes, localized at each longitudinal location x. GAWBS (Guided Wave Acoustic Brillouin 
Scattering)^"!! is a special case of this, in the low-frequency limit. Since neither Raman 
nor Brillouin excitations are completely localized, this treatment requires a frequency and 
wave-number cut-off, so that the field operator \1/ is slowly varying on the phonon scattering 
distance scale. The corresponding coupled nonlinear operator equations are: 



v lTx + m)^ x) 



—b(t,u,x) 



-Auj(x) 



u" d 2 



X E {x)&{t,x)$(t,x) 



2 dx 2 

R(u,x) b(t, to, x) + b\t, io, x) duo\^{t,x) 

L J 

iub(t, io, x) - iR(co, x)§ f (t, x)^(t, x) . 



(3.4) 



In summary, the original theory of nonlinear quantum field propagation is now extended 
to include both the the electronic and the Raman nonlinearities. The result is a modi- 
fied Heisenberg equation with a delayed nonlinear response to the field due to the Raman 
coupling. On integrating the Raman reservoirs, one obtains: 
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where 



-Au(x) + 



uj" d 2 



2 dx 2 Jo- 



dt' x (f, x) [ (t-t',x) + r (t, x) 



x (t, x) = x (x)S(t) + 26(t) / R 2 (u>, x) s\n(ut)duj 

Jo 

r R (t,a;) = - / R(u,x) b(t,u,x) + b\t,uj,x) dw, 
Jo 1 J 



(3.5) 



in which we have defined Q(t) as the step function. 

The operators T R and F R ' are stochastic, with Fourier transforms defined using the 
normal Fourier transform conventions for field operators: 



f R (u,x) 
f R \u,x) 



2tt 
1 



dtexp(iujt)f R (t,x) 
dtexp(-iut)T R (t,x). 



(3.6) 
(3.7) 



The frequency-space correlations are given by: 

(f *V, x >) f R (u, x)) = 2 x "(x, \u\) [n th (\cu\) + Q(-cu)} 5(x - x')5(cu - u'). (3.8) 

In this expression, we introduce a Raman amplitude gain of x" P er photon flux, equal 
to the imaginary part of the Fourier transform of x(t,x), so that: x"{ x -> 1^1) = kR 2 (x, \u>\). 
Here we use the Bloembergen normalization for response function Fourier transforms, 



x(uj,x) = J dtexp(iujt)x(t,x) , 



(3.9) 



which does not have the \/2ir factor included. 

Of some significance is the physical interpretation of the correlation functions, which can 
be regarded as directly contributing to the normally ordered spectrum of the transmitted 
field. Given a cw carrier, the correlations when u is positive correspond to an anti-Stokes 
(blue-shifted) spectral term, which is clearly zero unless the thermal phonon population is 
appreciable. However, when uj is negative, the theta function term indicates that the Stokes 
(red-shifted) spectral term is nonzero, due to spontaneous Stokes photons emitted even at 
zero temperature. 



B. Raman gain measurements 

The measured intensity gain due to Raman effects at a given relative frequency u per unit 
length, per unit carrier photon flux Jo = v is: 

1 9 In 7 ..... , . 

-2 X "(uj,x)/v 2 . (3.10) 



I dx 



Here the gain is positive for Stokes-shifted frequencies (u < ), and negative for anti-Stokes 
{uj > 0), as one would expect. This relationship allows the coupling to be estimated from 
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measured Raman gain and fluorescence properties. The simplest way to achieve this goal is 
to expand the Raman response function in terms of a multiple-Lorentzian model, which can 
then be fitted to observed Raman fluorescence data using a nonlinear least-squares fit. We 
therefore expand: 

n 

X (t,x) = x E {x)5{t) + x(x)6(t) Y l Fj8 j exp(-5 j t) sin(^t) . (3.11) 

3=0 

For normalization purposes, we have introduced x{ x )i which is defined as the total effective 
nonlinear phase-shift coefficient per unit time and photon density (in units of rad.m/s), 
obtained from the low-frequency nonlinear refractive index. This is given in terms of the 
electronic or fast-responding nonlinear coefficient, x E ( x ) > together with the Raman contri- 
bution, by integrating over time: 

/•oo roo 

x {x) = x E {x) + 2 / R 2 (uj,x)sm(ut)dujdt. (3.12) 
Jo Jo 

In the above expansion, Fj are a set of dimensionless Lorentzian strengths, and uij and 
Sj are the resonant frequencies and widths respectively, of the effective Raman resonances 
at each frequency. To improve convergence, the Lorentzian strength parameters are not 
constrained to be positive. The j = Lorentzian models the Brillouin contribution to the 
response function. In general, all of these parameters could be x-dependent, but we will 
often assume that they are constant in space for notational simplicity. The values for an 
n = 10 fit in the case of a typical fused silica fiber are given in Table [l], including an estimate 
of the effective Brillouin contribution averaged over the individual Brillouin scattering peaks. 
The coefficient of the electronic nonlinearity is now obtained explicitly in terms of the total 
nonlinear refractive index: 

X E (*) = ft(1 "^ W ° V . (3.13) 

where Uq is the carrier frequency, A is the effective cross-sectional area of the traveling mode, 
and / is the fraction of the nonlinearity due to the Raman gain, which has been estimated 
using the procedure outlined above: 



•yR 2 r 00 f°° 
f = — = — / dt dujR 2 (u;, x) sin(otf) 
X X J o Jo 

-0.2. (3.14) 

A result of this model is that the phonon operators do not have white noise behaviour. 
In fact, this colored noise property is significant enough to invalidate the usual Markovian 
and rotating-wave approximations, which are therefore not used in the phonon equations. 
Of course, the photon modes may also be in a thermal state of some type. However, thermal 
effects are typically much more important at the low frequencies that characterize Raman 
and Brillouin scattering, than they are at optical frequencies. In addition, if the input is 
a photon field generated by a laser, any departures from coherent statistics will be rather 
specific to the laser type, instead of having the generic properties of thermal fields. 

Finally, there is another effect which has been so far neglected. This is the ultra-low 
frequency tunneling due to lattice defects!!. As this is not strictly linear, it can not be 



S 



included accurately in our macroscopic Hamiltonian. Despite this, the effects of this 1/f 
type noise may be included approximately for any predetermined temperature. This can 
be achieved by simply modifying the refractive index perturbation term so that it becomes 
Au(t,x) and generates the known low-frequency refractive index fluctuations. 



4. Gain and absorption 

In silica optical fibers, there is a relatively flat absorption profile, with a minimum absorption 
coefficient of approximately 0.2dB /km in the vicinity of the commonly used communications 
wavelengths of around A = 1.5/im. This effect can be compensated for by the use of fiber laser 
amplifiers, resulting in nearly zero net absorption over a total link that includes both normal 
and amplified fiber segments. In practical terms, this situation leads to an approximately 
uniform fiber environment, provided the net gain and loss are spatially modulated more 
rapidly than than the pulse dispersion length. These additional effects need to be included 
within the present Hamiltonian model, in order to have a fully consistent quantum theory. 
For wide-band communications systems, either with time-domain multiplexing or frequency- 
domain multiplexing, it can become necessary to include the frequency-dependence and 
spatial variation of the gain and loss terms. This is especially true if spectral filters are 
included in the fiber. 



A. Absorbing reservoirs 

The absorption reservoir is modeled most simply by a coupling to a continuum of harmonic 
oscillators at resonant frequency uj. In the interaction picture used here, the Hamiltonian 
term causing rapidly varying operator evolution of the reservoir at the carrier frequency u 
is subtracted, leaving: 

POO fOC f ^ *l 

H' A = hj dx J du)\^(x)a i (u},x)A(u,x) + h.c. ] + {uj — u )(Fa(u,x)j , (4.1) 

where A(uj,x) provides the frequency dependent coupling between the radiation modes and 
the absorption reservoirs. The reservoir annihilation and creation operators, a and a) , have 
the commutation relations 

a(uj,x),a*(uj',x') = S(x - x')5(uj - uj 1 ) . (4.2) 

The equations for the absorbing photon reservoirs can be integrated immediately. The 
photon reservoir variable, for instance, obeys: 

8 ^ 

— a(t, cu, x) = —i(uj — uo)a(t, cu, x) — iA(u, x)^(t, x) . (4.3) 
Hence, the solutions are: 

a(t,u,x) = a(to,u, x) exp[— i(u — uo)(t — 1 )] — iA{u, x) \ exp[— i(u — u)o){t — t')\${t', x)dt' , 

J t 

(4.4) 
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with initial correlations for the reservoir variables in the far past (£ — > — oo) given by: 
{a* (to, u, x)a(to, uj', x')) = n t f l (u)S(x — x')8(u — oo') , 

(a(t , uj, x)a\t , uj', x')) = {n th (uj) + l]5(x — x')5(uj — uj') . (4.5) 

The solution for a(t, uj, x) is substituted into the Heisenberg equation for the field evolu- 
tion, giving rise to an extra time- dependent term, of the form: 

rOO POO ft ^ 

-i A*(u,x)a(t,uj,x)dw = - / dw\A(u, x)\ 2 / dt' exp[-i(u - u )(t - t')]^(t', x) 

JO JO Jt 

rco 

— ij dwA*(uj, x) exp[— i(u — uj )(t — t )]a(t , uj, x) 
Jo 

roo ^ ^ 

= - dt" 1 A (t",x)t!(t-t",x) + T A (t,x), (4.6) 
Jo 

where t" — t — t' , and the response function and reservoir terms are obtained most simply by 
extending the lower limit on the frequency integral to — oo, introducing only an infinitesimal 
error in the process, so that: 

/+oo 
duo\A(uj,x)\ 2 exp[-i(cu - u )t] 
-oo 

T A (t,x) = — i / du)A\u, x) exp[— i(uj — uj )(t — t )]a(to, uj, x). (4.7) 
Jo 

The response function integral represents a deterministic or 'drift' term to the motion, 
with a Fourier transform of: 

*I A (uj,x) = J 'j A (t,x)exp(iut)dt = 7 A (cj,x) + i^ A " '(to , x) , (4.8) 
so that the amplitude loss rate is: 

1 A (uj,x) = n\A(u + uj,x)\ 2 . (4.9) 

In the case of a spatially uniform reservoir with a flat spectral density, the Wigner-Weiskopff 
approximation (neglecting frequency shifts) gives a uniform Markovian loss term with: 

7 A (t)^7 A 5(t) , (4.10) 

where the average amplitude loss coefficient is: 

roo /"+oo 

y 4 = y^o) = / / dtduj \A(uj)\ 2 exp[-i{uj - u )t} 

Jo J -oo 

= 7 A + ^7 A ". (4.11) 

This approximation, known as the Markov approximation, is generally rather accurate for 
the absorbing reservoirs, whose response does not typically vary fast with frequency. An 
exception to this rule would be any case involving resonant impurities in the fiber, or very 
short pulses whose bandwidth is comparable to the frequency-scale of absorption changes. 
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The second quantity in Eq. ( |4.7|) , T A (t,x), behaves like a stochastic term due to the 
random initial conditions. Neglecting the frequency dependence of the thermal photon 
number, the corresponding correlation functions are 

^ ^ poo 

(T A (t,x)T A \t',x')) = / duj\A(uj,x)\ 2 exp\-i(u-u )(t-f)][n th (u) + l]5(x-x') 

Jo 

« [ 7 A (t - t', x) + 7 A *(t' - t, x)][n th (cu ) + 1] S(x - x') , (4.12) 

and: 

^ ^ poo 

(T A ^(t',x')T A (t,x)) = / duj\A(uj,x)\ 2 exp{-i(uj-uj )(t-t')}n th (Lu)5(x-x') 

Jo 

« [ 7 A (t - x) + ! M {t' - t, x)]n th {u ) S(x - a/) . (4.13) 

At optical or infra-red frequencies, it is a good approximation to set n t h(uJo) = 0. On Fourier 
transforming the noise sources, one then obtains: 

(t A (u, x)T a \uj\ x')) = 2 7 A (w, x)5(x - x')8(u - J). (4.14) 

Again taking the simplifying case of a spatially uniform reservoir in the Wigner-Weiskopff 
limit, this reduces to: 

(f A (t, x)f A t(t', x')) = 2 7 A 5(t - t')5(x - x') 

(T A \t,x)t A (t',x')) =0. (4.15) 

Note that the dimensions for the amplitude relaxation rates are [7 ] = s -1 . It is easy 
to show that 2^ A /v corresponds to the usual linear absorption coefficient for fibers during 
propagation. A typical measured absorption figure in current fused silica communications 
fibers is 0.2dB/km in the minimum region of absorption (near A = 1.5/im). The corre- 
sponding absorption coefficient is 2 r y A /v ~ 2.3 x 10~ 5 m _1 . To the extent that this effect 
is wavelength (and hence frequency) dependent, the resulting dispersion can be included as 
well, giving rise to a complete response function 7 A (t) for absorption. Non-Markovian ef- 
fects like this can either be neglected completely - which is a good approximation for slowly 
varying absorption in undoped fiber - or else included in the correlation functions of the 
reservoirs as given above. 

The physical meaning of the reservoir operator spectral correlations is best understood by 
considering the effect of these terms on photodetection, which according to photodetection 
theory, means a normally ordered field correlation. This involves normally ordered reservoir 
correlations to lowest order. Since these are zero, we conclude that the absorbing reservoirs 
essentially add no quantum noise that is observable via normal photodetection methods. 

B. Waveguide laser amplifiers 

The equations for gain or laser reservoirs are generally more complex, involving the nonlinear 
response of atomic impurities added to provide some gain in the fiber medium. This also 
involves a pump process (usually from a semiconductor laser) to maintain the lasing atoms 
in an inverted state. In the case of silica fibers, a commonly used lasing transition is provided 
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by erbium impurities^. The effect of these gain reservoirs is typically to introduce new types 
of dispersion, owing to the frequency dependence of the gained. In addition, there are new 
nonlinear effects, due to the effects of saturation - which in turn depend on the pumping 
intensity. 

It is possible to develop a detailed theory of erbium laser amplifiers. However, this 
paper will treat the simplest possible quantum theory of a traveling-wave quantum-limited 
laser amplifier. More details of the quantum theory, including nonlinear effects, are treated 
elsewhere^. However, the simple theory presented here provides a microscopic justification 
for observed quantum noise effects in fiber amplifier chains. In particular, it reproduces 
the results of the phenomenological theory of Gordon and HausB, which is known to give 
predictions in accord with soliton transmission experiments. The resulting "Gordon-Haus 
jitter" can be reduced through the use of filtering techniques. Assuming that the laser 
amplifier is polarization-insensitive, we again omit the polarization index. The reservoir 
variable cf^ = 1 1)^(2^ is an atomic transition operator, which induces a near-resonant atomic 
transition from an upper to a lower state, with two-level transitions having an assumed 
density of p(u, x) in position and resonant angular frequency to . 

These quantum effects are modeled here by including gain reservoirs in the Hamilto- 
nian, coupled by a frequency dependent term G(u,x) to the radiating field. Here the gain 
terms a ± (u;, x, t) represent the raising and lowering Pauli field operators, for two-level lasing 
transitions at frequency u. In more detail, we have gain given by an interaction Hamiltonian: 

H' G = h J°° dxJ^du^[m + {uj,x)G{uj,x) + h.c.} + U} ~ 2 LJ0 (r z (u,x) >l j , (4.16) 

where the atomic raising and lowering field operators, a ± , are defined in terms of discrete 
Pauli operators, by: 



IH ' 



a + (uj,x,t) = ^2\2)(l\ ll e}cp(-iu t)8(x - x^5{u - u t 

Sjp(uj,x) p 

a~(u,x,t) = 1 ^\l)(2\ fl exp(iu t)5(x - x^)5(uj - uj^) , 
sjp(u,x) p 

a z (u,x,t) = £[|2><2| " |1><1|U(* " ~ ^) • (4-17) 

These operators are in general time-dependent in the Heisenberg picture, but have the 
equal-time commutation relations: 

d- + (t,u,x),a~(t,u',x') = a z ( y t,cu,x)5(x - x')8(u - u'). (4.18) 

In the limit of complete inversion, with linear response and pure inhomogeneous broadening, 

d _ 

— a (t, u, x) = — i(w — ojQ)d (t, u>, x) + ia z (t, u>, x)G(u, x)^/(t, x) . (4.19) 

C/ v 

Hence, the solutions in the amplifier case are: 
<j~(t, u, x) = cr~(t , u, x) exp[— i(u — uj ){t — to)] 

+ iG{u,x) f exp[-i(u - u ){t -t')]a z (t'u,x)§(t',x)dt' . (4.20) 

Jto 
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With complete inversion, (d z (to, uj, x)) = 1 , so the initial correlations for the reservoir 
variables in the far past (to — > — oo) are given by: 

(cr + (t , uj, x)a~(t , to', x')) = S(x — x')5(uj — uj') , 

(a~ (t , uj, x)a + (t , uj', x')) = . (4.21) 

We substitute the solution for a~(t,uj,x) into the Heisenberg equation for the field evo- 
lution, assuming no depletion of the inversion, and trace over the atomic gain reservoirs. 
This gives rise to an extra time-dependent term, of the form: 



roo roo rt ^ 

-i G*(uj,x)a~(t,uj,x)duj = / duj\G(uj, x)\ 2 / dt' exp[-i(uj - uj )(t - t')]^(t', x) 

JO JO Jt 

roo 

— i dujG*(uj, x) exp[—i(uj — uj )(t — t )]a~ (t ,uj, x) 
Jo 

roo 

= dt" 1 G (t",x)^(t-t",x) + T G (t,x), (4.22) 
Jo 

where t" — t — t' , as before. This gives: 



/+oo 
duj \G(uj,x)\ 2 exp[-i(uj - uj )t] 
-oo 

^ rOO 

T G (t,x)^-i dujG*(uj,x)exp[-i(uj - uj )(t - t )]cr~(to,oj,x) . (4.23) 

J — oo 

Fourier transforming the response function gives: 

7 G (cj, x) = J 7 G (t, x) exp(iut)dt = >y G (uj, x) + ii G "(uJ, x) , (4.24) 
and the (real) resonant amplitude gain coefficient is: 

1 G (uj,x)=ii\G(uj + uJo,x)\ 2 . (4.25) 

As with the loss case, T G (t,x) behaves like a stochastic term due to the random initial 
conditions. The corresponding correlation functions are 

^ ^ roo 

(r G Ut',x')T G (t,x)) = / duj\G(uj,x)\ 2 exp[i(uj-uJo)(t-t')}5(x-x > ) 

Jo 

= [ 7 G (t - t', x) + 7 G * (t' - t, x)} S(x - x) . (4.26) 
Fourier transforming these noise sources gives: 

(f G V, x')f G (uj, x)) = 2 7 G (cu, x)5(x - x')5(uj - uj') (4.27) 

Taking the uniform fiber in the Wigner-Weiskopff limit as before, so 7 G = 7 G (0,x), this 
reduces to: 

(f G (t,x)f G ^(t',x')) =0 

{f G ^t,x)T G (t',x')} = 2 1 G 5(t-t')5(x-x'). (4.28) 
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The dimensions for the amplitude gain are [y G } = s _1 . On Fourier transforming, the 
response function can be related to the measured intensity gain 2Re[7 G (c<;, x)/v ] at any 
frequency offset u, relative to the carrier frequency uj . This allows one to obtain the linear 
gain coefficient for fibers during propagation. Since measured laser gain figures can be much 
greater than the absorption, it is possible to compensate for fiber absorption with relatively 
short regions of gain. The results presented here are only valid in the linear gain regime. 
More generally, a functional Taylor expansion up to at least third order in the field would be 
needed to represent the full nonlinear response of the laser amplifier, together with additional 
quantum noise terms. 

Finally, it is necessary to consider the result of incomplete inversion of an amplifier. 
Here, the noninverted atoms give rise to absorption, not gain, and will generate additional 
quantum-noise and absorption response terms. These must be treated as in the previous 
section, including non-Markovian effects if the absorption line is narrow-band. An important 
consequence is that the measured gain only gives the difference between the gain and the 
loss. This doesn't cause any problems with the deterministic response function - but it 
does cause difficulties in determining the amplifier quantum noise levels, which can only be 
uniquely determined through spontaneous fluorescence measurements. Obviously, the lowest 
quantum noise levels occur when all the lasing transitions are completely inverted. 

The physical meaning of the reservoir operator spectral correlations for the amplifier 
case is clearly quite different to the case of the absorber. If we consider the effect of these 
terms on photodetection as before, which means a normally ordered field correlation, we 
should look again at the normally ordered correlations of the reservoirs. Since these are 
no longer zero, we conclude that the amplifying reservoirs emit fluorescent photons due to 
spontaneous emission over the amplifier bandwidth. 



5. Combined Heisenberg Equations 



Coupling linear gain and absorption reservoirs in this way to the Raman-modified Heisenberg 
equation leads to a generalized quantum nonlinear Schrodinger equation. Such equations are 
sometimes called quantum Langevin equations. In the present case of a single polarization, 
the resulting field equations are: 

/ d d \ - r°° 

[v— + — J *(t, x) = -J o dt' 7 (t', x)*(t - t', x) + F(t, x) 

i 

dt' x (t')[m}(t-t',x) + r R (t,x) 



+ 1 



uj" d 2 

2 dx 2 Jo 



In this equation, 

7 (t, x ) = 7 A (t, x) - 7 G (t, x) + iAiu(x)5(t) 



$(t,x). (5.1) 



(5.2) 



is a net linear response function due to a coupling to linear gain/absorption reservoirs, 
including the effects of a spatially varying refractive index. This can be Fourier transformed, 
giving: 7(0;, x) = j(u,x) + ij"(u>,x), where 7(0;, x) < for gain, and 7(0;, x) > for 
absorption. Similarly, T(t,x) is the linear quantum noise due to gain and absorption. The 
actual measured intensity gain at frequency u> + u is given in units of [m _1 ] , by: 
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^ = 2( 7 G (^, x) - 7 A (^, x))/v . (5.3) 
The stochastic terms have the correlations 

(f *V, x') Y R (u, x)) = 2 X "(x, \u\)[n th (\cj\) + e{-u)]6{x - x')8(u - u') 
(f V, x') f(w, x)> = 2 7 g (cj, x)8(x - x')8(u - w') 

(T{uj, x) f V, x')) = 2 7 A (w, x)8{x - x')8{uj - </) , (5.4) 

where we have introduced minimal linear quantum noise terms Y and for the 
gain/absorption reservoirs, and where thermal photons have been neglected (since usually 
Tiojq » kT, as explained in section |j). Equation (|5.1|) can be easily generalized to include 
nonlinear absorption or laser saturation effects, relevant to amplifiers with intense fields, but 
these terms are omitted here for simplicity. 

This complete Heisenberg equation gives a consistent quantum theoretical description 
of dispersion, nonlinear refractive index, Raman/GAWBS scattering, linear gain, and ab- 
sorption. It is important to notice that the reservoir correlations have a simple physical 
interpretation, especially in the zero-temperature limit. Normally ordered noise correlations 
occur when there is gain, anti-normally ordered correlations when there is absorption. This 
is the reason why the normally ordered Raman noise correlations vanish at zero tempera- 
ture for positive frequencies. At low temperatures, Raman processes only cause absorption 
to occur at positive detunings from a pump frequency. Thermal correlations have a more 
classical behaviour, and occur for both types of operator ordering. 

It is often useful to do calculations in a simpler model, in which we include the effects of 
uniform gain and loss in a moving frame. This can either be carried out using a standard 
moving frame (x v = x — vt), or with a propagative time (t v = t — x/v ) as in the original 
Gordon-Haus calculations. For propagative calculations, it is most convenient to use photon 
flux operators 

$(t„,x) = y/v$(t,x). (5.5) 

For long pulses, assuming a uniform gain/loss response in the frequency domain, the prop- 
agative transformation gives the following approximate equations: 



^-$(t v , x) = - r < ^^$(4 - c x) + m/v^ 

ox Jo v 



+ i 



2 du 2+ Jo- v 2 [ l[v v ' } + v 



${t v ,x). (5.6) 



In addition, if the pulses are narrow-band compared to the gain and loss bandwidths, and the 
reservoirs are uniform, then the gain and absorption reservoirs are nearly delta-correlated, 
with 

(ft(t, x v ) Y(t', <)) = 2 1 G 8(x v - x v ')5(t - t') 

(Y(t, x v ) Y\t\ x' v )) = 2 1 A 8{x v - x v ')8(t - t'). (5.7) 

It is essentially this set of approximate equations that corresponds to those used to 
predict the solitonEil self-frequency shiftH and related effects@ in soliton propagation, except 
for the omission of the Raman reservoir terms. 
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6. Phase-space methods 



The Heisenberg equations are not readily soluble in their present form. To generate nu- 
merical equations for analytic calculations or for simulation, operator representation the- 
ory can be used. There is more than one possible method, depending on which phase- 
space representation is used. The positive-P representation, for example, produces exact 
provided phase-space boundary terms are negligible, while a truncated Wigner 
represent ationS^ gives approximate results that are valid in the limit of large photon num- 
ber. It is important to note that the Wigner method represents symmetrically ordered rather 
than normally ordered operator products, and so has finite quantum noise terms even for 
a vacuum field. These can be thought of as corresponding to the shot noise detected in a 
homodyne or local-oscillator measurement, while the positive-P representation represents 
normally ordered operators, and therefore corresponds to direct-detection noise. 

Either technique can be used for this problem, each with its characteristic advantages and 
disadvantages. The positive-P representation, although exact, uses an enlarged phase-space 
which therefore takes longer to simulate numerically. It only includes normally ordered noise 
and initial conditions, and this is an advantage in some cases, since the resulting noise is 
zero in the vacuum state. The Wigner technique is simpler, and for large mode occupations, 
its results are accurate enough for many purposes. However, it has the drawback that it 
includes symmetrically ordered vacuum fluctuations. 

First, we expand the field operators in terms of operators for the free-field modes. Apply- 
ing the appropriate operator correspondences to the master equation for the reduced density 
operator in which the reservoir modes have been traced over, namely 

p 9 = TT R p = Tr R — \H } p] , (6.1) 
in 1 J 

gives a functional equation for the corresponding operator representation. 

In the positive-P case, the equation is defined on a functional phase-space of double the 
classical dimensions, so that a complete expansion in terms of a coherent-state basis is 
obtained: 

fa® = J J p(t,y,v) l -^dmdm. (6.2) 

The resulting Fokker-Planck equation for the positive distribution P(t, \I>, \1>) has only second 
order derivative terms. Sometimes the notation = \1/ is used to indicate the stochastic 
field that corresponds to the hermitian conjugate of 

The equation for the Wigner function W(t, \P) also contains third and fourth order deriva- 
tive terms, which may be neglected at large photon number. The resultant Fokker-Planck 
equation in either case, can be converted into equivalent Ito stochastic equations for the 
phase space variables \I> (and ^). Physical quantities can be calculated by forming the 
average over many stochastic realizations, or paths, in phase-space. For example, in the 
positive-P representation, (\& ^) stochastic = (^vp) quantum, while in the Wigner representa- 
tion, (^*#) stochastic = + ^^quantum- 
It should be clear from this that the positive-P representation directly generates an 
intensity corresponding to the usual normally ordered intensity that is detected in direct 
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photodetection. The Wigner representation, however, generates an intensity result that 
includes some vacuum fluctuations. In a computer simulation with a finite number M of 
modes, we must correct the Wigner result by subtracting M/2 from any simulated photon 
number, or vM/2 from any calculated photon flux, in order to obtain the direct photodetec- 
tion result. For the calculation of a homodyne measurement, the Wigner method will give 
the most directly suitable result with symmetric ordering. In this case it is the positive-P 
representation that will need correction terms added to it. Once these corrections are made, 
either method will give similar results, although the sampling error may not be identical. 



A. Modified nonlinear Schrodinger equation 



Standard custom in fiber optics applications^ involves using the propagative reference frame 
with the normalized variables: r = (t — x/v)/t and ( = x/x , where t is a typical pulse 
duration used for scaling purposes, and x = tl/\k"\ ~ 1km for dispersion shifted fiber. This 
change of variables is useful only when slowly varying second order derivatives involving ( 
can be neglected, which occurs for vto/xo <C 1. For typical values of the parameters used 
in experiments, this inequality is often well-satisfied (vto ~ 10 _4 m). To make it simpler to 
compare with this usage, we will make the same transformation for the stochastic equations 
that are equivalent to our complete operator equations, and scale the variables used in a 
dimensionless form. 

For definiteness, we will now focus on the spatially uniform case. The resultant equation, 
which includes gain and loss, is a Raman-modified nonlinear Schrodinger (NLS) equation 
with stochastic noise terms: 



d 



(t,C) 



dT's(T-TW,C) + r(T,o 

dr '/,( r _ r ')^( r ' )C )0( r ' )C) + r«(r,C) 



i d 2 (f) f°° 



(r,0, 



(6.3) 



where = ^yvt /n is a dimensionless photon field amplitude. The photon flux is \cf)\ 2 n/t , 
and n = \k"\Ac/ (n2hu 2 to) = v 2 t /x%o is the typical number of photons in a soliton pulse 
of width to? for scaling purposes. The positive sign in front of the second derivative term 
applies for anomalous dispersion ik" < 0), which occurs for longer wavelengths, and the 
negative sign applies for normal dispersion (k" > 0). A similar equation is obtained in the 
positive-P case, except that </>*and T R *(r, () are replaced by non-complex-conjugate fields 
denoted (fi + and T r+ (t, () respectively: 



C) = - / dr'g*(r - r')0 + (r', fl + r + (r, C) 
i d 2 ( 



+ 



2 dr 2 



- i / dr'h*(r - tW, 0<f> + (r', + 



+ (r,C), (6.4) 



The equations in and <p + both have the same additive noises and identical mean values, 
only differing by the independent parts of the multiplicative noise sources - which therefore 
generate nonclassical quantum statistics. 
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The causal linear response function g{r) is defined as: 
g[r) = . (6.5) 

V 

If the Fourier transform of this function is g(Q) = g(fl) + ig'(Q), then we can relate this to 
dimensionless intensity gain a G (Q) and loss a A (Q), at a relative (dimensionless) detuning 
of fl , by: 

2g(tt) = a A (Q) - a G (Q) . (6.6) 

The causal nonlinear response function h(r) is normalized so that / h{r)dr = 1, and it 
includes both electronic and Raman nonlinearities: 

h(r) = h E (r) + h R (r) = ^^o) (6J) 

The Raman response function h R {r) causes effects like the soliton self-frequency shiftH. The 
response function Fourier transform is given by: 

h(n) = J dtexp(inr)h(T) = ti{ti) +ih"(tt). (6.8) 

This definition has the property that the value of h(Q) = h(uito) is a dimensionless number, 
which depends on the frequency u only, independent of the time-scale used for normalization. 
The Raman gain, whose spectrum has been extensively measured^!, can be modeled as a 
sum of n Lorentzians, as explained in section |3] and as illustrated in Fig. 1. 
This expansion as n Lorentzians gives a response function of the form 

n 

h R (t/t ) = e(t)^F i 5 i t exp( —6jt)sm((jjjt), (6.9) 

3=0 

It is most convenient to express these in terms of dimensionless parameters flj = ujjto and 
Aj = 6jt , giving: 

n 

h R (r) = 0(r) FjAj exp(-A i r) sin(fi J r). (6.10) 

j=o 

Here Aj are the equivalent dimensionless widths (corresponding to damping), and the flj 
are the dimensionless center frequencies, all in normalized units. It is useful to compare 
these results with the dimensionless Raman gain a R (Q) , normalized following Gordon!!, 
which uses a characteristic time-scale of to . The relationship of macroscopic coupling R(ou) 
to measured Raman gain a R (fl) is R 2 (u) = x& R {uto)/2ii. It follows that the dimensionless 
gain function is 

a R (n) = 2\h"(Q)\ . (6.11) 

These stochastic partial differential equations can be discretized and, without any fur- 
ther approximation, can be numerically simulated!!!! using a split-step Fourier integration 
routine. The equations include all the currently known noise physics significant in soliton 
propagation, including effects like the soliton self-frequency shift. Guided acoustic wave 
Brillouin scatteringil'Bi! noise sources are included in the Raman gain function. These 
have little effect on the position of an isolated soliton, but are important for long-range 
soliton collision effects!! that occur in pulse-trains. 
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B. Initial conditions 



The initial conditions for the calculations could involve any required quantum state, if the 
+P representation is used. In the case of the Wigner equations, only a subset of possible 
states can be represented with a positive probability distribution. The usual initial condition 
is the multi-mode coherent state, since this is the simplest model for the output of mode- 
locked lasers. In general, there could be extra technical noise. We note that the choice of 
a coherent state is the simplest known model of a laser sources. To represent this in the 
positive-P distribution is simple; one just takes 

P (r,O) = [0£(r,O)]* = (0(r,O)>. (6.12) 

In the Wigner case, which corresponds to symmetric operator ordering, one must also include 
complex quantum vacuum fluctuations, in order to correctly represent operator fields. For 
coherent inputs, the Wigner vacuum fluctuations are Gaussian, and are correlated as 

(<Mr,o)) = (0M)> 

(A<Mt, O)A0^(t', 0)) = ^5(r - t'). (6.13) 

in 

We note that these equations imply that an appropriate correction is made for losses at the 
input interface, so that the mean- field evolution is known at the fiber input face. 



C. Wigner noise 

Both fiber loss and the presence of a gain medium contribute quantum noise to the equations 
in this symmetrically ordered representation. The complex gain/absorption noise enters the 
Wigner equation through an additive stochastic term T, whose correlations are obtained 
by averaging the normally and anti-normally ordered reservoir correlation functions given 
previously, together with appropriate variable changes. This symmetrically ordered noise 
source is present for both gain and loss reservoirs. Thus, 

(r(n, c)r*(ir, c')> = 1 1 J 2 I_ 1 )} 6(c - C)5(n - of), (q.u) 

where T(fl, () is the Fourier transform of the noise source: 
1 r°° 

T(Q,C) = -= dTF(r, C) exp(iQr) 

V J-oo 

r*(fi,C) = 4= drr(r,C)exp(-ifir). (6.15) 

V ITC J-oo 



Similarly, the real Raman noise, which appears as a multiplicative stochastic variable T R , 
has correlations 

(r R (n, c)r**(Q', O) = ^-P^ \n th (\n\/t ) + 1} s(( - (')S(n - n') . (e.ie) 
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Thus the Raman noise is strongly temperature-dependent, but it also contains a sponta- 
neous component which provides vacuum fluctuations even at T = 0. Since the spontaneous 
component can occur through coupling to either a gain or a loss reservoir, in a symmetrically 
ordered representation, it is present for both positive and negative frequency detunings. 

It must be remembered here that the noise terms in the Wigner representation do not 
correspond to normally ordered correlations, and so have no direct interpretation in terms 
of photodetection experiments. Any predictions made with this method of calculation need 
to be corrected by subtracting the appropriate commutators, to convert the results into a 
normally ordered form. This is the reason why there is no obvious distinction between the 
amplifier and absorber cases. 



D. -\-P noise 



The positive P-representation is a useful alternative strategy, because it does not require 
truncation of higher order derivatives in a Fokker-Planck equation, and corresponds directly 
to observable normally ordered, time-ordered operator correlations. It has no vacuum fluctu- 
ation terms. Provided the phase-space boundary terms are negligible, one can then obtain a 
set of c-number stochastic differential equations in a phase-space of double the usual classical 
dimensions. These are very similar to the classical equations. Here the additive stochastic 
term is as before, except it only depends on the gain term a G ; the conjugate term T* is used 
in the + equation: 

(r(n, C) r*(fi', c')> = -^<KC - 0*(« - n'). (6.17) 

n 

Since this representation is normally ordered, the only noise sources present are due to 
the gain reservoirs. There is no vacuum noise term for the absorbing reservoirs, because 
absorption simply maps a coherent state into another coherent state. 

The complex terms T R , T R+ include both Raman and electronic terms (through h'(£l)). 
As elsewhere in this paper, we regard r iJ+ (r2, Q as a hermitian conjugate Fourier transform 
(with the opposite sign frequency exponent): 

r + (fi,C) = -= / rfrr + (r,C)exp(-^r). (6.18) 

V ATX J—oo 



. 2tt 

This quantity is not the same as r fi *(f2, C), since it involves a noise source that is in general 
independent. In some cases, where classical noise is dominant (and nonclassical squeezing 
is negligible), we can ignore this fact, and approximately set T R+ (£l,() — T i? *(fi,C)- More 
generally, we obtain the following results: 

(T R (Q, C) r*(Q', O) = *(C - CO S(Q + ST) {[n th (\n\/t ) + 1/2] a R (\n\) - i h'(Q)} In 

(r R+ (n', CO r*(n, 0) = 5(( - CO 5(n - n') [n th (\n\/t ) + e(-n)] a R (\n\)/n. (6.19) 

This equation is an expected result, since it states that when Q < the spectral intensity 
of noise due to the Stokes process, in which a photon is down-shifted in frequency by an 
amount fl with the production of a phonon of the same frequency, is proportional to n th + 1. 
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However the anti-Stokes process in which a phonon is absorbed (fi > 0) is only proportional 
to n th . Thus, at low temperatures the only direct noise effect is that due to the Stokes 
process, which can be interpreted physically as originating in spontaneous photon emission, 
detectable through photodetection. 

As one might expect, the two forms of equation are identical at high phonon occupation 
numbers, when classical noise is so large that it obscures the differences due to the operator 
orderings of the two representations. Another, less obvious, result is that the two equations 
have identical additive noise sources, provided the gain and loss are balanced. To understand 
this, we can see that in the absence of any net gain or loss, the differences in the operator 
correlations due to ordering is a constant, contained in the initial conditions. 

However, when gain and loss are not equal, the additive noise sources are quite different. 
In particular, the Wigner representation has noise contributions from both types of reservoir. 
On the other hand, the normally ordered +P method only leads to additive noise when 
there is a real fluorescent field present, which is detectable through photodetection. This 
corresponds physically to some kind of gain, either due to the presence of an amplifier, or 
through Raman effects. 

In general, the Wigner and +P reservoir correlations are obtainable simply by examining 
the expectation values of the Heisenberg reservoir terms, with symmetric and normal order- 
ing respectively. The additional term proportional to h'(Q) in the +P noise terms is due 
to dispersive nonlinear effects, and gives rise to a nonclassical noise source which is respon- 
sible for the observed quantum squeezing effects. Extensions required to treat polarization 
dependent Raman scattering are given elsewhereB. 



7. Conclusions 

Our major conclusion is that quantum noise effects due to the intrinsic finite-temperature 
phonon reservoirs and finite bandwidth amplification or absorption can be readily modeled 
using stochastic equations. The equations themselves have the usual classical form, together 
with correction terms that we can describe as quantum noise terms. The precise form of the 
correction terms depends in detail on the representation employed (although this difference 
is purely due to operator ordering), as well as the physical origin of the reservoir couplings. 
These correction terms can be non-Markovian or nonuniform in space. The generation of 
the corresponding stochastic noises is a straightforward numerical procedure, and generally 
much simpler than the use of noncommuting operators. By contrast, the original operator 
equations have no practical numerical solution in most cases, due to the exponential growth 
of the dimension of the underlying Hilbert space with the number of modes and photons 
involved. 

Detailed applications to short-pulse soliton communications will be given in a subsequent 
paper!. In general, the increasing bandwidth, reducing pulse energies and greater demands 
placed on fiber communications and sensors mean that these quantum limits are becoming 
increasingly important. Already, limits set by quantum amplifiers are known to have great 
significance in long-distance laser-amplified communications systems. We note that the 
quantum theory given here also establishes the levels of quantum noise in silica fibers in more 
general situations. Examples of this are for dispersion-managed fiber communicationsS'0, 
and for fiber ring lasers with relatively low gainta. Similarly, these equations set the limits 
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for experiments using spectral filtering and related techniques to generate sub-shot- 
pulsesHS in optical fibers. 
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FIGURES 

Fig. 1. The parallel polarization Raman gain |9{/i(u;to)}| = \h"(uto)\ for the 11-Lorentzian model 
(continuous lines) and the single-Lorentzian model (dashed lines), for a temperature of T = 300K. 
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TABLES 

Table 1. Fitting parameters for the 11-Lorentzian model of the Raman gain function h R (t/to). All 
frequencies are in T.rad/s. 



7 
J 


3 


fi); 

3 




o 


0.16 


0.005 


0.005 


l 


-0.3545 


0.3341 


8.0078 


2 


1.2874 


26.1129 


46.6540 


3 


-1.4763 


32.7138 


33.0592 


4 


1.0422 


40.4917 


30.2293 


5 


-0.4520 


45.4704 


23.6997 


6 


0.1623 


93.0111 


2.1382 


7 


1.3446 


99.1746 


26.7883 


8 


-0.8401 


100.274 


13.8984 


9 


-0.5613 


114.6250 


33.9373 


10 


0.0906 


151.4672 


8.3649 



27 



